Bose glass transition and spin-wave localization for 2D bosons in a random potential 
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A spin-wave (SW) approach of the zero temperature superfluid — insulator transition for two 
dimensional hard-core bosons in a random potential /i = ±W is developed. While at the classical 
level there is no intervening phase between the Bose-condensed superfluid (SF) and the gapped 
disordered insulator, the introduction of quantum fluctuations lead to a much richer physics. Upon 
increasing the disorder strength VK, the Bose-condensed fraction disappears first, before the SF. 
Then a gapless Bose-glass (BG) phase emerges over a finite region, until the insulator appears. 
Furthermore, in the strongly disordered SF regime, a mobility edge in the SW excitation spectrum 
is found at a finite frequency Qc, decreasing with VK, and presumably vanishing in the BG phase. 



A correct understanding of the interplay between 
strong correlations and disorder is one of the most diffi- 
cult question in condensed matter physics [1, 2]. While 
Anderson theory of localization [3] for single particle 
states is now a well-established paradigm to describe elec- 
tronic transport in disordered environments, the equiv- 
alent bosonic problem of dirty superconductors or su- 
perfluids remains quite challenging [4, 5]. Despite nu- 
merous pionneer studies [4, 6], several questions remain 
open. For instance in ID the universal character of the 
Luttinger exponent at the SF-BG transition has been 
recently questionned [7-9]. For more realistic higher 
dimensional systems, relevant for disordered supercon- 
ductors [10, 11], quantum antiferromagnets [12], or cold 
atoms [13], quantum Monte Carlo approaches have con- 
siderably improved our understanding of the dirty boson 
problem all along the past twenty years [14-19], but have 
also raised new issues regarding the universal value of 
some critical exponents [18-21], and so far have only ad- 
dressed ground-state properties. On the analytical side, 
important progresses have been made recently to go be- 
yond mean- field (MF) theory [22-25]. Although a naive 
MF is unable to find a localization transition, even at 
very strong disorder [6], a quantum cavity approach on 
the Bethe [22] or the square lattice [23, 25, 26] is able 
to capture such a transition. Nevertheless, several is- 
sues remain unsolved, in particular concerning finite fre- 
quency physics [22, 27, 28], and the outstanding question 
of many-body localization [29-31]. 

In this letter, we want to improve our understanding of 
the interplay between quantum fluctuations and disorder 
by addressing the spin- wave (SW) corrections for the Ma- 
Lee model in a disordered bimodal potential, defined on 
the square lattice by the Hamiltonian 



^b = -tE(%+^i^O 



(y) 



W y^ eiUi 



(1) 



which describes preformed Cooper pairs (hard-core 
bosons) hopping between nearest neighbor sites with a 
random chemical potential VKe^, where e^ = ±1 with 
probability 1/2. In the disorder- free case (i.e. e^ = 1 
for instance), this well-known model [32, 33] displays 



two phases at T = : (i) a Bose-condensed super- 
fluid regime for incommensurate filling < (n) < 1 if 
\W\ < 4t, and (ii) a trivial insulator, filled (empty) with 
(n) = 1 ((n) = 0) for W > U {W < -4t). Using the 
Matsubara-Matsuda mapping [34] of hard-core bosons 
onto pseudo-spin 1/2, Hamiltonian (1) is exactly equiv- 
alent to a spin-^ XY model in a longitudinal field along 
the z axis. A mean-field description, where spin oper- 
ators are treated as classical vectors with two angles Oi 
and (/)i, gives an energy £ = — | ^^^^-^ sin ^^ sin Oj cos(^i — 

= constant and 



.) 



2 ^{ij) 

fyjj ^ ^- Ci COS Oi^ minimized by 
cos^^ = eiW/{4:t) if W < 4t, meaning XY order for the 
spins (and superfluid Bose condensate for the bosons). If 
W > 4:t there is no XY order anymore: all spins point 
along the z axis with cos^^ = e^ which, in the bosonic 
language, corresponds to a disordered insulator with lo- 
cal occupations (n^) = (1 + ei)/2 (= or 1). In the XY 
regime, condensate and superfluid densities (po and psf) 
are both equal to (sin^ Oi)/4 = [l - (W/4t)^] /4, vanish- 
ing dit W = At. Within such a classical description, a 
direct transition between SF and gapped phases is ob- 
served for W > 4t, as visible in Fig. 1, with no interme- 
diate localized regime, an artifact of MF theory. 

However, when quantum fluctuations are introduced, 
the situation changes dramatically [35]. Before describ- 
ing in more details our SW results, let us first briefly 
summarize our main conclusions. Here, we have studied 
square systems up to 64 x 64 for several hundreds of 
disordered samples, which allowed us to get infinite size 
extrapolations for various thermodynamic quantities 
such as the superfluid and the condensate densities 
Psf and pQ. An intervening gapless Bose glass phase 
is unambiguously found between the superfluid and 
the gapped insulator. Properties of the SW excitation 
spectrum have also been studied, namely, the sound 
velocity and the inverse participation ratio (IPR) of 
the SW excited states [36-38]. The localization of SW 
modes displays very interesting features vs frequency Q. 
We find a finite mobility edge f^c, such that states with 
frequencies Q < Qc are extended and states at Q > Qc 
are localized. Upon increasing the disorder strength, ftc 
decreases and presumably vanishes in the BG phase. 
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FIG. 1: (color online) Superfluid (SF) and Bose-Einstein con- 
densate densities psi and po plotted together with the gap A, 
against the disorder strength W/t. The classical densities (•) 
both vanish at the same point W — At whereas SW corrected 
quantities pf^ (D) and p^ (o) vanish at different points 
Wo < VKsf < 4t, leaving a finite window for an intervening 
gapless Bose glass before the gapped insulator. Insets depicts 
SF and insulating phases in the pseudo-spin representation. 
Disorder average was performed over several hundreds of dis- 
ordered samples. The green line is a guide to the eyes. 



Let us now present in more details these results. SW 
corrections for hard-core bosons are treated in a straight- 
forward way [32, 33], first making a local rotation for the 
pseudo-spin operators, and then introducing Holstein- 
Primakoff bosons (a,a'*'). At the linear SW level, the 
hard-core bosons model (1) reads Hb = S ^1-L^'^\ with 



^(2)^ 



where t,- 



--^[(ii^a-aj- 
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n, - t[e,ej{v/Atf - 1], 
with V = max(iy, 4t) and u = min(lV,4t). Because 
translationnal invariance is broken by the disorder, the 
quadratic bosonic Hamiltonian Eq. (2) is diagonalized by 
a Bogoliubov transformation in real space which yields 



N 



1. 



H(^)=^b,(ata, + -) 



(3) 



ftp are the SW frequencies and {a ^a^) describe Bogoli- 
ubov quasi-particles. In the clean case W = 0, the modes 
p are labeled by the wave vectors k = ikx^ky) and the SW 
spectrum fl\^ = It^J A — 2 (cos k^ + cos ky) ~ 2t|k| when 
|k| -^ 0, recovering the linear Bogoliubov spectrum with 
a "velocity of sound" vq = 2t. 

It is important to note that Bose-condensate and 
superfluid fractions are intrinsiqually different objects 
which are only equal in the simplest ME description; ^He 
being one of the best examples of a strongly correlated 



(non ME) bosonic system with po/Psf — 8% at low tem- 
perature [39]. To go beyond ME, we want to compute the 
first SW corrections for the condensate and the superfluid 
response. As discused in detail in Ref. 33, there are two 
ways for correctly computing 1/S corrections to a physi- 
cal observable O. One may evaluate its expectation value 
(O) in the l/5'-corrected ground-state, but this is not an 
easy task for our disordered problem. Perhaps more sim- 
ply one can add a small symmetry-breaking term to the 
Hamiltonian of the form SH = — FO, compute the 1/S- 
corrected energy and take the derivative with respect to 
the field in the limit F ^ 0. For instance, the condensate 
density po = -^^ijiK^j)^ which is simply related in 
the pseudo-spin language to the transverse magnetization 



(Po 



when N ^ oo) is obtained by adding a term 



— F ^- Sf to the pseudo-spin XY model. The SF density 
Psf can be equally computed using the response of the sys- 
tem to twisted boundary conditions [40] , via the helicity 
modulus (or superfluid stiffness) Tgf = d'^E{(f)/d(f'^\ ^q, 
then simply related to the SF density by psf = ^sff^t. 

Numerical results for (S^) on lattices up to 32 x 32 
(averaged over several hundreds of disordered samples), 
are shown in panel (b) of Fig. 2 versus 1/N. There, we 
clearly see that when the disorder exceeds W/t = 3.5, 
SW correction starts to become larger than the classi- 
cal contribution, thus giving a negative magnetization 
which we interpret as a transition to a zero magnetization 
state. Finite size extrapolations to the thermodynamic 
limit [full lines in Fig. 2(b)] give the disorder average con- 
densate density po = ((^^))^ [41], plotted in Fig. 1. Such 
a behavior is not surprising as it is well-known that quan- 
tum fluctuations on top of the classical solution deplete 
the condensate mode. Here quantum fluctuations coop- 
erates with disorder, leading to a monotonous destruction 
of Bose condensation, gradually increasing from ~ 25% 
depletion at W = up to 100% at Wo/t = 3.55(5). 
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FIG. 2: (color online) (a) Crossing of the disorder average su- 
perfluid density psf x L^ in the vicinity of the critical disorder 
where superfluidity disappears. Using z = 2, a very con- 
vincing crossing is found for Wsi = 3.738(1). (b) Transverse 
magnetization (S^), computed using a numerical derivative 
with respect to a small transverse field F = t/100 and aver- 
aged over several hundreds of samples, plotted against 1/A^ 
for various disorder strengths in the vicinity of Wo = 3.55(5). 



More surprising is the behavior of the SF density psf , 
computed in the presence of a smah twist angle Lp = 10~^. 
Infinite size extrapolations for psf are shown in Fig. 1 
(blue squares) where we see that contrary to the conden- 
sate, quantum fiuctuations first enhance superfiuidity for 
weak disorder, until W/t = 3 where quantum and dis- 
order effects start to cooperate and destroy the super- 
fiuid which finally disappears for a critical disorder 4 > 
Ws^/t = 3.738(1) > Wo/t. Interestingly, one can also test 
hyperscaling at the 2D critical point where [15] Psf ^ L~^ 
is expected. As shown in Fig. 2 (a) we check that the best 
crossing of psf x L^ is obtained at Wsf/t = 3.738(1) with 
a critical exponent z = 2.0(1), in a surprisingly good 
agreement with the expected z = d [4, 18]. 

Before discussing the properties of the excitations, let 
us briefiy comment the fact that Bose condensation and 
superfiuidity disappear for different values of the disor- 
der: Wq < VFsf- While this result appears potentially 
quite interesting, a careful caution is required to interpret 
it. Indeed, the fact that for Wq < W < Wsi quantum 
fiuctuations have destroyed the Bose condensed fraction 
whereas a finite SF density remains, would mean that a 
condensate-free superfiuid has been stabilized, a concept 
recently discussed in the context of the frustrated anti- 
ferromagnet BaCuSi206 at large fields [42]. Nevertheless 
here, contrary to the frustration-induced case where the 
absence of condensate is quite fragile since it relies on a 
special symmetry of the lattice, it is legitimate to wonder 
whether the window Wgf — ^o remains finite beyond lin- 
ear SW corrections, a question perfectly suited to future 
quantum Monte Carlo simulations [43]. In any case, we 
have demonstrated here that linear SW corrections can 
drive a bosonic state where both po and psf are zero over 
a finite window W G [ VFgf , 4t] , which is interpreted as a 
insulating Bose glass with a gapless excitation spectrum, 
as we discuss now. 

We first focus on the first excitation level above the Bo- 
goliubov vaccum. We find the entire regime < W/t < 4 
to be gapless, with a zero mode Qq c^ 0, and a finite size 
gap to the first excitated state scaling in the limit L ^ 1, 
as Asw(I/) ^ 27rv/L^ as visible in Fig. 3 (a) for various 
values of the disorder W. The prefactor v is identified 
with the velocity of sound (or SW velocity) and is shown 
in Fig. 3 (b), rescaled by its zero-disorder value vq = 2t, 
versus W/t. In the same panel (b), the classical hydro- 
dynamic relation for the velocity v = Y^Tgf//^, is also 
plotted, Tgf and hz being the MF results for the helicity 
modulus and the compressibility. Both estimates for v 
compare remarkably well. Interestingly, the bottom of 
the SW spectrum is only weakly affected by the disor- 
der and remains phonon-like (delocalized) over the entire 
gapless regime W/t G [0,4] with a finite velocity, almost 
disorder-independant, except very close to the insulating 
phase at W/t = 4 where v abruptely drops down [44]. 
This finite velocity in the entire gapless regime is con- 
sistent with recent studies of Anderson localization of 
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FIG. 3: (color online) (a) Finite size SW gap plotted vs 1/L 
for various disorder strengths VK < 4t in the gapless regime. 
Full lines are quadratic fits of the form A(L) = 27vv/L-\-b/L'^, 
where v, the sound velocity, is displayed in panel (b) against 
W/t, together with the estimate from the classical hydrody- 
namic relation (see text). The full blue line is a power-law fit 
- (4t - W)^-^^^. (c) Infinite size extrapolation of the SW gap 
(red circles) in the insulating regime W > 4t. The full red 
line is a power- law fit ^ (W — 4t)°"^^^, and the black dotted 
line is the clean case (fi = W > 4t) result: A = VK — 4t. 



phonons in disordered solids [45, 46]. Above W = At the 
zero mode disappears and a finite gap opens in the SW 
spectrum, as visible in Figs. 1 and 3 (c). Interestingly, 
this gap does not scale linearly with 1^— 4t as in the clean 
case, but opens up more rapidly, presumably ~ y/W — 4t 
and approaches the clean case only at large W. 

Following [37], we have investigated the localization 
properties of the entire SW Bogoliubov excitation spec- 
trum. Here we shall just mention the main results of this 
study which will be described in details in a longer arti- 
cle [38]. In Ref. 37, it has been observed that the localiza- 
tion properties of the SW excited states depend crucially 
on the frequency in a way similar to the Anderson local- 
ization of phonons [45]. Here, we have analyzed this effect 
by considering the inverse participation ratio (IPR), de- 
fined for each (normalized) state \p) = X^^^^l^), where 
i are lattice sites, by IPRp = J^i^i l<^?l^- ^^^ delocal- 
ized modes IPR ^ 1/N whereas localized states display 
a finite IPR ~ 1/^^, where ^ is the localization length. 
Since SW spectra are discrete for finite size systems, in 
particular at low energy, we define disorder average IPRs 
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FIG. 4: (color online) Inverse participation ratio IPR in the strongly disordered SF phase for W = 3 At; (a) Best crossing of 
IPRx Ar^2/2 obtained with N = 256, • • • , 1600 at a mobility edge Qc/vo ^ 1.2 with a fractal dimension D2 = 1.48. (b) IPRx AT 
plotted vs N for different frequencies Qi = r^/'^o (different symbols). The dashed red line ^ N shows the fully localized case 
when N ^ ^^ and the full black line is the critical scaling ^ jy^-^2/'^ ^^^ ^j^g mobility edge. 



over finite slices of frequencies centered around Q: 



iFR{n) 



ZG{np,n±SQ) 



(4) 



where e{np,n±6n) = lifft-Sn <np< n^Sn, andO 
otherwise, with SQ/vq = 1/20 in the following. While for 
weak disorder W/t < 2 all the excited states are found de- 
localized, similarly to the clean case where the coefficients 
are simply the Fourier modes a^ = exp(ikp • Yi)/\/N ^ 
thus giving for all frequencies IPR(r^) x TV = (9(1), the 
case of strongly disordered SF appears much more in- 
teresting, as visible in Fig. 4 which shows representa- 
tive results for W/t = 3.4. At low energy the modes 
are delocalized, but the situation changes dramatically 
above a certain threshold frequency Vtc - the mobility 
edge - where IPR(r^) x N starts to increase linearly 
with A^, a characteristic signature of localization. At 
the mobility edge, as in the case of the Anderson tran- 
sition [1, 45], the IPR is found to display an anoma- 
lous scaling IPR(17c) oc A/'~^2/2 ^ith a fractal dimension 
D2 — 1.48 < 2. This is well visible in the panel (a) of 
Fig. 4 where the best crossing of IPR((^) x A/'^2/2 ^^^ 
been obtained for D2 = 1.48. Note that for other disor- 
der strengths, the same fractal exponent has been found 
to obtain the best crossing curves, separating extended 
modes ai Q < Qc from localized ones ai Q > Qc- 

The evolution of the mobility edge Qc against increas- 
ing disorder is shown in Fig. 5 where we see that l^c ^ 
when the BG phase is approached. While the localization 
transition point is easily identified in Fig. 4 for W/t = 3.4, 
closer to the SF-BG boundary error bars for Qc get big- 
ger. Indeed, it becomes more difficult to correctly esti- 
mate the localization transition on finite size systems for 
W/t > 3.7 since the crossing point displays a sizeable 
drift towards smaller frequencies when N increases. This 
observation supports the fact that the BG phase is local- 
ized for all r^ > 0. The phase diagram energy - disorder 
in Fig. 5 displays 3 different regimes: (i) delocalized exci- 



tations in the SF regime below a finite mobility edge Qc] 
(ii) absence of modes below a finite gap A for W/t > 4; 
(iii) localized excitations above Qc or A. Finally one can 
mention that, contrary to Refs. 22, 47, inside the insulat- 
ing phase, we do not find any mobility edge from local- 
ized excited states at small frequency to extended states 
at large frequencies. Conversely, our results support the 
idea that superfluidity emerges out of the localized BG 
phase by a delocalization at (^ = 0, in agreement with 
Refs. 27, 37 

It remains to consider to what extent these results de- 
pend on the particular distribution of the disorder used 
and the SW approximation made. Also, larger system 
sizes should allow to draw firm conclusions regarding the 
absence of a mobility edge inside the Bose-glass phase. 
Work is in progress along such lines [38, 43]. 
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FIG. 5: (color online) Energy of the SW excitations Q (in 
units of vo) as a function of W/t. All states are extended 
(delocalized) below the mobility edge Qc and localized above. 
The shaded area represents the localized-delocalized bound- 
ary with quite large error bars close to the SF-BG transition 
point where we expect Qc -^ 0. In the gapped insulating side, 
there is no state below the gap A, and all excitations above 
are localized, and connected to other localized excitations. 
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